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ABSTRACT 

We investigate the linear tidal perturbation of a viscous Keplerian disk by a companion star orbiting 
in a plane inclined to the disk. We consider m = 1 perturbations with odd symmetry with respect to 
the z = midplane. The response frequency may be either finite or vanishing. These long-wavelength 
perturbations produce a well-defined warp. 

Since the response of a viscous disk is not in phase with the perturbing potential, a tidal torque is 
exerted on the disk. When the perturber rotates outside the disk, this torque results in a decrease of the 
disk angular momentum, and thus in an increase of its accretion rate. 

We show that this tidal torque is comparable to the horizontal viscous stress acting on the background 
flow when the perturbed velocities in the disk are on the order of the sound speed. If these velocities 
remain subsonic, the tidal torque can exceed the horizontal viscous stress only if the viscous parameter 
a v which couples to the vertical shear is larger than the parameter ah coupled to the horizontal shear. 
In protostellar disks, bending waves, which are predominantly excited in the outer regions, are found to 
propagate and transport a significant fraction of the negative angular momentum they carry deep into 
the disk inner parts. 

If the waves are reflected at the center, resonances occur when the frequency of the tidal waves is 
equal to that of some free normal global bending mode of the disk. If such resonances exist, tidal 
interactions may then be important even when the binary separation is large. Out of resonance, the 
torque associated with the secular perturbation, which is proportional to a v , is generally much larger 
than that associated with the finite frequency perturbations. As long as the waves are damped before 
they reach the center, the torque associated with the finite frequency perturbations does not depend on 
the viscosity, in agreement with theoretical expectation. These calculations are relevant to disks around 
young stars and maybe also to disks in X-ray binary systems. 

Subject headings: accretion, accretion disks — binaries: general — hydrodynamics — 
stars: pre-main sequence — waves 



1. INTRODUCTION 

1.1. Tidal Interactions in Pre-main- Sequence Binary 
Systems: Coplanar versus Noncoplanar Systems 

Thanks to recent optical images taken with the Hub- 
ble Space Telescope, accretion disks around low-mass 
pre-main sequence stars known as T Tauri are no 
longer objects of theoretical speculations (McCaughrean & 
O'Dell 1996). As expected, they appear to be common 
around young stars (Stauffer et al. 1994). Furthermore, 
most T Tauri stars are observed to be in multiple systems 
(see Mathieu 1994 and references therein). Since the dis- 
tribution of the separation of pre-main sequence binaries 
has a peak around 30 astronomical units, which is smaller 
than the extent of the disks in such systems (Edwards et 
al. 1987), it is expected that many circumstellar disks will 
be subject to tidal effects due to the influence of binary 
companions. 

The theory of tidal interaction has already success- 
fully satisfied many observational tests. For example, it 
has been predicted that angular momentum exchange be- 
tween orbital motion and disk rotation through tidal in- 
teraction may truncate the disk (see, e.g., Papaloizou & 



Pringle 1977; Paczynski 1977). This has recently received 
observational confirmation, both through the study of 
young binaries spectra (Dutrey, Guilloteau & Simon 1994; 
Jensen, Mathieu & Fuller 1996) and direct imaging (Rod- 
dier et al. 1996). 

The tidal effect of an orbiting body on a differentially 
rotating disk has been well studied in the context of plane- 
tary rings (Goldrcich & Tremaine 1978), planet formation, 
and interacting binary stars (see Lin & Papaloizou 1993 
and references therein). In these studies, the disk and orbit 
are usually taken to be coplanar. However, there are obser- 
vational indications that the plane of the disk and that of 
the orbit may not necessarily be aligned. The most strik- 
ing evidence for such noncoplanarity is given by HST and 
adaptive optics images of the pre-main sequence binary 
HK Tau (Stapelfeldt et al. 1998). The projected binary 
separation in this system is 340 AU, the observed radius 
of the circumsecondary disk is 105 AU and the disk as- 
pect ratio at a radius of 50 AU is about 0.08. The authors 
suggest that the primary would be located at an angle of 
at least 20° above the disk plane. In the case of the hier- 
archical triple system Ty Cra, some models suggest that 
the orbit of the close spectroscopic binary and that of the 
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wider binary are not in the same plane (see Corporon, La- 
grange & Beust 1996; Bibo, The & Dawanas 1992). Also, 
observations of the pre-main sequence binary T Tau have 
revealed that two bipolar outflows of very different orien- 
tations originate from this system (Bohm & Solf 1994). 
Since it is unlikely that they are both ejected by the same 
star, each of them is more probably associated with a dif- 
ferent member of the binary. Furthermore, jets are usually 
thought to be ejected perpendicularly to the disk from 
which they originate. These observations then suggest 
that disks and orbit are not coplanar in this binary sys- 
tem. More generally, observations of molecular outflows in 
star forming regions commonly show several jets of differ- 
ent orientations emanating from unresolved regions of ex- 
tent as small as a few hundred astronomical units (Davis, 
Mundt & Eisloffcl 1994). 

We note that the spectral energy distribution of a cir- 
cumstellar disk could be significantly affected by the lack 
of coplanarity in binary or multiple systems (Terquem & 
Bertout 1993, 1996). This is because reprocessing of radi- 
ation from the central star by a disk depends crucially on 
the disk geometry. 

In such noncoplanar systems, both density and bending 
waves are excited in the disk by the perturbing companion. 
They are respectively of even and odd symmetry with re- 
spect to reflection in the disk mid-plane. Since much work 
has already been carried out on density waves, which arc 
the only waves excited in coplanar systems, we shall focus 
here on bending waves. Because their wavelength is larger 
than that of density waves, they are expected to propagate 
and transport angular momentum further radially into the 
disk. 

1.2. Bending Waves 

Bending waves were first studied by Hunter & 
Toomrc (1969) in the context of the Galactic disk. They 
calculated the dispersion relation associated with these 
waves taking into account self-gravity but ignoring pres- 
sure and viscosity. Bending waves were further con- 
sidered in the context of galactic disks (Toomrc 1983; 
Sparke 1984; Sparke & Casertano 1988) and planetary 
rings (Shu, Cuzzi & Lissauer 1983; Shu 1984). In these 
studies, the horizontal motions associated with the ver- 
tical displacement of the disk were neglected. However, 
Papaloizou & Pringlc (1983) showed that these motions 
are important in gaseous Keplcrian disks, where they are 
nearly resonant. These authors pointed out that when the 
disk is sufficiently viscous so that its response to verti- 
cal displacements is diffusive rather than wave-like, these 
motions lead to the decay timescale of the warp being 
times the viscous timescale t„ ;V - Here a v is the Shakura & 
Sunyaev (1973) parameter that couples to the vertical 
shear, and t v , y is the disk viscous timescale associated 
to a v . When a v is smaller than unity, this timescale is 
smaller than t vy , which is the decay timescale of the warp 
obtained by previous authors. 

In the context of protostellar disks, Artymowicz (1994) 
studied the orbital evolution of low-mass bodies crossing 
disks, and Ostriker (1994) calculated the angular momen- 
tum exchange between an inviscid disk and a perturber on 
a near-parabolic noncoplanar orbit outside the disk (see 
also Heller 1995 for numerical simulations). The vertical 



structure of the disk was not taken into account in these 
studies. Papaloizou & Lin (1995) performed a full 3 dimen- 
sional linear analysis of m = 1 bending waves in accretion 
disks. Comparing analytical results with numerical simu- 
lations, they showed that a vertical averaging approxima- 
tion can be used when the wavelength of the perturbation 
is much larger than the disk semi-thickness. Using the 
WKB approximation, they calculated the dispersion re- 
lation associated with bending waves taking into account 
self-gravity, pressure and viscosity (see also Masset & Tag- 
ger 1996 for a study of thick disks). They found that in 
a non self-gravitating thin Keplerian disk in which a v is 
small compared to the disk aspect ratio, bending waves 
propagate without dispersion at a speed which is about 
half the sound speed. Using these results, Papaloizou & 
Terquem (1995) calculated the m = 1 bending wave re- 
sponse of an inviscid disk with the rotation axis misaligned 
with a binary companion's orbital rotation axis. This 
work was extended by Larwood et al. (1996) who per- 
formed nonlinear calculations of the disk response using 
smoothed particle hydrodynamics (SPH). This numerical 
method was further used by Larwood & Papaloizou (1997) 
who analyzed the response of a circumbinary disk subject 
to the tidal forcing of a binary with a fixed noncoplanar 
circular orbit, and by Larwood (1997) who considered the 
tidal interaction occurring between a disk and a perturber 
following a parabolic trajectory inclined with respect to 
the disk. 

The present work should be contrasted with earlier 
studies of purely viscous disks in which pressure and 
self-gravity were ignored. The behavior of nonplanar 
purely viscous accretion disk subject to external torques 
has been considered previously by a number of authors. 
Following the analysis of Bardeen & Petterson (1975), 
these studies have been developed in connection with 
the precessing disks of Her X-l and SS 433 for instance 
(e.g., Petterson 1977; Katz 1980a, 1980b), galactic warps 
(e.g., Steiman-Cameron & Durisen 1988), or warped disks 
around AGN and in X-ray binaries (Pringle 1996). How- 
ever, the response of a purely viscous disk is very different 
from that of a viscous gaseous disk, which is a better repre- 
sentation of a protostellar disk. In the former case, evolu- 
tion occurs only through viscous diffusion, whereas in the 
latter case pressure effects manifesting themselves through 
bending waves can control the disk evolution. When these 
waves propagate (i.e. when they are not damped), they 
do so on a timescale comparable to the sound crossing 
time. Since this is much shorter than the viscous diffu- 
sion timescale, communication through the disk occurs as 
a result of the propagation of these waves. Even when the 
waves are damped before reaching the disk boundary, a dif- 
fusion coefficient for warps that is much larger than that 
produced by the kinematic viscosity may occur because of 
pressure effects (Papaloizou & Pringle 1983). Self-gravity 
may also strongly modify the response of viscous disks, 
because like pressure it allows bending waves to propagate 
and its effect is global. For a more detailed discussion of 
the respective effects of viscosity, pressure and self-gravity, 
see Papaloizou, Terquem & Lin (1998). 

1.3. Torque Exerted on the Disk by a Companion on an 
Inclined Orbit 
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There can be no tidal torque exerted on an inviscid disk 
with reflective boundaries that contains no corotation reso- 
nance (Goldreich & Nicholson 1989). However, if the disk 
is viscous or has a nonreflective boundary, the compan- 
ion exerts a torque on the disk that leads to an exchange 
of angular momentum between the disk rotation and the 
orbital motion. We can distinguish the following three 
components of the torque: 

1. The component along the line of nodes m akes the 
disk precess at a rate that can be uniform (see § |5.2 and 



Kuijken 1991; Papaloizou & Terquem 1995; Papaloizou et 
al. 1998). Although we will calculate in this paper the pre- 
cession frequency induced by the tidal torque, we will not 
discuss the disk precession in detail since this has already 
been done in the papers cited above. 

2. The component of the torque along the axis per- 
pendicular to the line of nodes in the disk plane in- 
duces the evolution of the inclination angle of the disk 
plane with respect to that of the orbit. This evolution, 
which is not necessarily toward coplanarity, occurs on the 
disk's viscous timescale if the warp is linear (Papaloizou & 
Terquem 1995; Papaloizou et al. 1998). Since it affects the 
geometry of the system only on such a long timescale, we 
will not discuss this effect here. We note that both these 
components of the tidal torque modify the direction of the 
disk angular momentum vector. 

3. The component of the torque along the disk's spin 
axis does not modify the direction but the modulus of the 
disk's angular momentum vector, thus changing the disk's 
accretion rate. This effect is important since it is a source 
of angular momentum transport in the disk, and it will be 
the main subject of this paper. 

So far, angular momentum exchange between a viscous 
gaseous disk rotation and the orbital motion of a compan- 
ion on an inclined circular orbit has not been considered. 
It is the goal of this paper to present an analysis of this 
process. In the case of an inviscid disk considered by Pa- 
paloizou & Terquem (1995), such angular momentum ex- 
change occurs because the waves are assumed to interact 
nonlinearly with the background flow before reaching the 
disk's inner boundary. Thus, only in-going waves propa- 
gate through the disk. This induces a lag between the disk 
response and the perturbation, which enables a torque to 
be exerted on the disk. In a viscous disk, the lag required 
for the torque to be non zero is naturally produced by the 
viscous shear stress. 



1.4. Plan of the Paper 

The plan of the paper is as follows. In § |2| we derive an 
expression for the perturbing potential. In § ^ we outline 
the basic equations. The model of the unperturbed disk 
is described in § |4[ The disk response is presented in § |[ 
where we give the linearized equations, study the disk re- 
sponse for both zero and finite perturbing frequencies, es- 
tablish the boundary conditions and describe the domain 
of validity of this analysis. In § |^ we calculate the torque 
exerted by the perturber on the disk, and relate it to the 
disk viscosity by using the conservation of angular mo- 
mentum. We show that this tidal torque is comparable to 
the horizontal viscous stress acting on the background flow 
when the perturbed velocities in the disk are on the order 
of the sound speed. If these velocities remain subsonic, the 



tidal torque can exceed the horizontal viscous stress only 
if the viscous parameter a v that couples to the vertical 
shear is larger than the parameter ah coupled to the hori- 
zontal shear. Numerical results are presented in § |?], both 
for uniform and nonuniform viscosity. They indicate that, 
for parameters typical of protostellar disks, bending waves 
are able to propagate and transport a significant fraction 
of the negative angular momentum they carry deep into 
the inner parts of the disk. They also show that, if the 
waves can be reflected at the disk's inner boundary, res- 
onances occur when the frequency of the tidal waves is 
equal to that of some free normal global bending mode 
of the disk. Out of resonance, if the disk and the or- 
bital planes are not too close to being perpendicular, the 
torque associated with the zero-frequency perturbing term 
is found to be much larger than that associated with the 
finite frequency terms. As long as the waves are damped 
before they reach the center, the torque associated with 
the finite frequency perturbations does not depend on the 
viscosity, in agreement with theoretical expectation (Gol- 
dreich & Tremaine 1982). We summarize and discuss our 
results in § |8[ 

2. PERTURBING POTENTIAL 

We consider a binary system in which the primary has 
a mass M p and the secondary has a mass M s . The binary 
orbit is circular with separation D. We suppose that the 
primary is surrounded by a disk of radius R <C D and 
of negligible mass M so that the orbital plane does not 
precess and the secondary describes a prograde Keplerian 
orbit with angular velocity: 



G (M a + M p ) 
D 3 



about the primary. In the absence of the secondary star, 
the disk is nearly Keplerian. We adopt a nonrotating 
Cartesian coordinate system (x, y, z) centered on the pri- 
mary star. The z-axis is chosen to be normal to the initial 
disk mid-plane. We shall also use the associated cylindri- 
cal polar coordinates (r, tp, z) . The orbit of the secondary 
star is in a plane that has an initial inclination angle 5 
with respect to the (x, y) plane. 

We adopt an orientation of coordinates and an origin of 
time such that the line of nodes coincides with, and the 
secondary is on, the cc-axis at time t — 0. We denote the 
position vector of the secondary star by D with D = |D|. 

The total gravitational potential \& due to the binary at 
a point with position vector r is given by: 



GM V 



D| 



GM s r ■ D 
D 3 ' 



(1) 



where G is the gravitational constant. The first dominant 
term is due to the primary, while the remaining perturbing 
terms are due to the secondary. The last indirect term ac- 
counts for the acceleration of the origin of the coordinate 
system. 

We are interested in disk warps that are excited by terms 
in the potential that are odd in z and have azimuthal mode 
number m = 1 when a Fourier analysis in (p is carried out. 
The terms in the expansion of the potential that are of the 
required form are given by: 
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(2) 



We consider a thin disk such that z <C r. Since r <C -D, 
we then retain only the lowest order terms inr/D and 2/£> 
in the above integral, which becomes: 



3 GM S 



4 D 3 



■rz [(1 — cos <5) sin 5 sin (<p + 2a>i) 



— ( 1 + cos <5) sin 6 sin (<p — 2wt) + sin 26 simp]. (3) 

Because the principle of linear superposition can be 
used, the general problem may be reduced to calculating 
the response to a complex potential of the form: 



\& = ifrze 



i(ip-Qpt) 



(4) 



Here the pattern frequency fl p of the perturber is one of 
0, 2uj or — 2o>, and / is an appropriate real amplitude. 



3. BASIC EQUATIONS 

The response of the disk is determined by the equation 
of continuity: 



dp 

dt 



+ V • (pv) = 0, 



and the equation of motion: 



9v 

dt 



+ (v V) v = - (-VP + F„) - V*, 



(5) 



(6) 



where P is the pressure, p is the mass density and v is 
the flow velocity. We allow for the possible presence of a 
viscous force but, following Papaloizou & Lin (1995), 
we shall assume that it does not affect the undistorted 
axisymmetric disk so that it operates on the perturbed 
flow only. This approximation is justified because the disk 
viscous timescale is much longer than the characteristic 
timescales of the processes we are interested in here (in 
particular, bending waves propagate with a velocity com- 
parable to the sound speed, so that the waves crossing time 
is much shorter than the viscous timescale). 

For simplicity, we adopt a polytropic equation of state: 



P = Kp 1+1/n , 



(7) 



where K is the polytropic constant and n is the polytropic 
index. The sound speed is then given by c 2 = dP/dp. 

4. EQUILIBRIUM STRUCTURE OF THE DISK 

We adopt the same model as in Papaloizou & 
Terqucm (1995). Namely we suppose that the equilibrium 
disk is axisymmetric and in a state of differential rotation 
such that, in cylindrical coordinates, the flow velocity is 
given by v = (0, rfl, 0) . For a barotropic equation of state, 
the angular velocity is a function of r alone. 

For a thin disk under the influence of a point mass, in- 
tegration of the vertical hydrostatic equilibrium equation 
gives (Papaloizou & Terquem 1995): 



HWH [2A ^ + """)' / "" t "(^") W1 ' 1 M 



p(r,z)=[2K(l + n)C%\ 



-n/(2n+l) 



/ 2 \ n 

x (m K ) 2 " /(2 " +1) (l - iL j . 



(9) 



Here is the Keplerian angular velocity given by £l 2 K = 
GMp/r 3 , H is the disk semi-thickness, S is the surface 
mass density defined by: 



and 



S(r) 



H(r) 



H(r) 



p (r, z) dz, 



r(i)r(n + i) 



where T denotes the gamma function. 

It then follows from the radial hydrostatic equilibrium 
equation that the angular velocity is given by: 



fi 2 = £l 2 K + 



2n 1 V(2«+l) 



2CZ 



d_ 

dr 



2/(2n+l) 



(10) 



In principle, the surface density profile should be de- 
termined self-consistently as a result of viscous diffusion 
(Lynden-Bell & Pringle 1974). However, since we ignore 
the effect of viscosity on the unperturbed disk, we are free 
to specify arbitrary density profiles for the equilibrium disk 
model. For S we then take a function that is constant in 
the main body of the disk with a taper to make it vanish 
at the outer boundary: 



S(r) = S - 



cr(r) 



which approximates £(r) = S min [1, cr(r)] with 

p> N 

o-(r) 



(11) 



1 - 



r- (R- A) 



The parameters A, which represents the width of the taper 
interior to the outer boundary, N, q and p are constrained 
such that the square of the epicyclic frequency, which is 
defined by: 

2 _ 20 d (r 2 Q) 
r dr 

is positive everywhere in the disk. For the numerical cal- 
culations we take N = (2n + 1) /2, A/R = 0.1, q = 4 and 
p=l. 

The polytropic constant K is determined by specifying 
a maximum value of the relative semi-thickness, H/r, in 
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the disk, and the fiducial surface mass density So is fixed 
by specifying the total mass M of the disk. 

Figure |l| shows the surface mass density, S, the rela- 
tive semi-thickness, H/r, the ratio of the sound speed to 
the Keplerian frequency, c s /CIk ~ H, and the ratio of 
the epicyclic frequency to the Keplerian frequency, k/CIk, 
versus the radius r for the equilibrium disk models de- 
scribed above with M = 10~ 2 M p and (H/r) max = 0.1. 
The epicyclic frequency departs from Qk at the outer edge 
of the disk because of the abrupt decrease of E there. 



5. THE DISK RESPONSE 
5.1. Linear Perturbations 

In this paper we consider small perturbations, so that 
the basic equations (|5|) and (0) can be linearized. As the 
relevant linearization of the equations has already been 
discussed by Papaloizou & Lin (1995), we provide only an 
abbreviated discussion here. We denote Eulerian pertur- 
bations with a prime. Because the perturbing potential, 
is proportional to exp [i (ip — dpi)] , the ip and t de- 
pendence of the induced perturbations is the same. We 
denote with a tilde the part of the Eulerian perturbations 
that depends only on r and z. For instance, the Eule- 
rian perturbed radial velocity is denoted v' r , and we have 
v' r (r, ip, z, t) = v' r (r, z)exp [i (p - Cl p t)]. 

We define a quantity g (the physical meaning of which 
will be given at the end of this section) through: 



pi 

h = -irzCl 2 g, 

P 



(12) 



and we also have g(r, ip, z, t) — g(r, z)exp [i (ip — Cl p t)]. 

It was first shown by Papaloizou & Pringle (1983) that 
horizontal motions induced in a near Keplerian tilted disk 
are resonantly driven. The main effect of a small viscos- 
ity is then to act on the vertical shear in the horizon- 
tal perturbed velocities dv' r /dz and dv' /dz. Demianski & 
Ivanov (1997) and Ivanov & Illarionov (1997) have recently 
given the relativistic generalization of this effect. Follow- 
ing Papaloizou & Lin (1995), we then retain only the hor- 
izontal components of the viscous force perturbation: 



'--U^y (13) 

where v v is the kinematic viscosity that couples to the ver- 
tical shear. This approximation is accurate when a v <C 
1. The numerical calculations have actually shown that, 
when a v becomes larger than 0.1, the terms that are ne- 
glected in this approximation become comparable to the 
terms that are retained. For this reason, we shall limit our 
numerical calculations to values of a v smaller than 0.1. 

Here we have assumed that the disk is turbulent and 
that dissipation can be modeled as a turbulent viscosity. 
The most likely mechanism for producing such a turbulent 
viscosity is the magnetorotational Balbus-Hawley instabil- 
ity (Balbus & Hawley 1991). We note however that this 
instability may not operate in all the parts of protostel- 
lar disks (Balbus & Hawley 1998). Also, it is not clear 
that the tidal waves would be damped by the turbulent 
viscosity in the way described by equations (fL3|). In this 



formulation, we implicitly assume that the turbulent vis- 
cosity is not affected by the tidally-produced velocities, 
which is not necessarily true. However, since the wave- 
length of the tidal waves we consider is much larger than 
the disk semi-thickness, this formulation may be correct. 

Papaloizou & Lin (1995) have shown that when the ra- 
dial wavelength of the perturbations is larger than the 
disk semi-thickness, a vertical averaging approximation, 
in which g is assumed to be independent of z, can be used. 
Each of v' r and v' v are then proportional to z, with v' z being 
independent of z. They have found some support for the 
procedure in a numerical study of the propagation of dis- 
turbances in disk models that take the vertical structure 
fully into account. 

Linearization of the equation of motion followed by mul- 
tiplication by pz and vertical averaging then gives the ve- 
locity perturbations in the form: 



it 

ciz 



(1-Cl p /Cl)g 



(i - cl p /cl) rf + g(3 - tyn) ci 2 /n 2 
(i-n p /n-ia v f - n 2 /n 2 



fj 



ciz i - ci p /ci V 1 2Ci 2 ciz 
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(14) 



(15) 



(16) 



rCl l-O p /0' 

where a v is a vertical average of the standard Shakura & 
Sunyaev (1973) parameter that couples to the vertical 
shear. It is defined by: 



pv v dz 



ci 



pz 2 dz 



(17) 



As mentioned above, the main effect of viscosity is to 
damp the resonantly driven horizontal velocities. There- 
fore, the viscous terms which do not act directly on the 
resonance have been neglected. This is consistent with 
keeping only the components of the viscous force given by 
equation (|13|). Again, this is accurate only when a v <C 1. 

The same procedure applied to the continuity equation 



gives: 



-i (ft -ftp) I- 



n 2 g 



(CI - Cl p ) 



(CI 



Or 



rfi- 



ci p y - Cl z 



i/j,- 



dr 



(18) 



in which 



J 



+°° pz 2 



dz, n = 



pz 2 dz 1 X = 



P ^'dz. 



We note that equation ( |18| ) depends on the viscosity only 
through the perturbed velocities. 

To close the system of equations, we need a relation 
between P' and p' . For simplicity, we assume that the 
equation of state is the same for the perturbed and un- 
perturbed flows. Linearization of the equation of state (|7]) 
then gives 
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P = c 2 J. (19) 

We now give a physical interpretation of g. Equa- 
tion ( |l6|) can be used to express g in terms of the vertical 
Lagrangian displacement £ z : 



2 £z 



ff = i(i-fip/n) - 



(20) 



where again we have defined £ 2 such that £ z (r, <p, z, t) = 
£ z (r,z)exp[i (ip - fl p t)}. 

When H p <C Ct, the physical relative vertical displacement 
is then given by: 



Re — R e (9) sm {<P — ^p*) + ^ TO (s) cos {<P ~ topt) ■ 

Thus, when fl p — 0, Re{g) and Im{g) are the relative ver- 
tical displacement along the y and x-axis respectively. We 
note that a constant value of g corresponds to a rigid tilt. 

5.2. Zero-Frequency Response 

From equation (|^), we have to consider the disk response 
to a secular potential perturbation with zero forcing fre- 
quency. 
(H) and 
differentia 



When ftp = and a v < 1, equations (|14|), (p.5|) 



reduce to the single second-order ordinary 
equation for g (Papaloizou & Lin 1995): 



d_ 

dr 



dg 



(l-iov) -k 2 /^ 2 rfr 
Here we have assumed that: 



r 



(22) 



where Slo = fi(-R) and r,„ is the disk inner boundary. Here 
the disk has been assumed to precess as a rigid body. This 
is expected if the disk can communicate with itself, ei- 
ther through wave propagation or viscous diffusion (self- 
gravity being ignored) , on a timescale less than the preces- 
sion period. As long as a v is smaller than the disk aspect 
ratio H/r (which is almost constant throughout most of 
the disk), bending waves may propagate (Papaloizou & 
Pringle 1983; Papaloizou & Lin 1995). The ability of 
these waves to propagate throughout the disk during a pre- 
cession time implies approximately that H/r > \ui p \ /O.Q. 
This is the condition for near rigid body precession in an 
inviscid disk (Papaloizou & Terquem 1995). The numer- 
ical simulations of Larwood et al. (1996) show that this 
condition also holds in a viscous disk in which a v does not 
exceed H/r. This indicates that in the absence of self- 
gravity, as long as a v < H/r (which is likely to be satisfied 
in protostellar disks), pressure is the factor that controls 
the efficiency of communication throughout the disk and 
then the precessional behavior of the disk. A more com- 
plete discussion of the precession of warped disks (which 
includes self-gravity) is given by Papaloizou et al. (1998). 

The Coriolis force can be taken into account by replac- 
ing with: 



\E f/ + 2iuj p rzQ,sin8. 



(25) 



H < Re < r— 
dr 



(23) 



(l-iov) -k 2 /Q 2 

which is certainly a reasonable approximation for a v « 1. 

It can be shown that with the term X in its present form, 
the x-component of the torque is non zero. This torque 
leads to the precession of the disk, as in a gyroscope. By 
writing the equation of motion for Q p — in the frame 
defined in § |^, we have assumed that the response of the 
disk appears steady in a nonrotating frame. Because of the 
precessional motion induced by the secular perturbation, 
the response is actually steady in a frame rotating with 
the precession frequency. The Coriolis force correspond- 
ing to this additional motion must then be added in the 
equation of motion (the centrifugal force can be neglected 
since it is of second order in the precession frequency). 
The magnitude of the Coriolis force in the frame in which 
the response appears steady has to be such that the x- 
component of the total torque (involving both the Coriolis 
and the gravitational forces) is zero. 

This requirement leads to the expression of the preces- 
sion frequency lo p (Papaloizou & Terquem 1995; see also 
Kuijken 1991 in the context of galactic disks): 



fio 




(24) 



For the zero-frequency case, ^' will represent the to- 
tal term (^5|) and the real amplitude / will represent 
/ + 2ujpfl sin 6 from now on. 

With this new definition of /, g is then given by: 



«(!■) / , - L "" ,> ' ,r K2/n2 I J[r")f[r")dr"dr', 

(26) 

where the quantity in the outer integral has to be evalu- 
ated at the radius r' . We have supposed that g = at 
r = ri ni but we note that an arbitrary constant may be 
added to g. This corresponds to an arbitrarily small rigid 
tilt that we assume eliminated by the choice of coordinate 
system. 

The above expression shows that, depending on whether 
a v is larger or smaller than 1 — n 2 /VI 2 ~ H 2 /r 2 , the imag- 
inary part of g is larger or smaller than its real part, re- 
spectively. According to equation ([jl|), this corresponds to 
a vertical displacement produced by the secular response 
being mainly along the x- or y-axis, respectively. 

We finally comment that the above treatment is accu- 
rate only when |wJ /fio is smaller than the maximum of 
K 2 /f2 2 ~ H 2 /r and a v . If this is not the case, then the 
precession frequency has to be taken into account in the 



1 



resonant denominator (1 
be the dominant term. 



ia v ) — k /CI where it would 
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5.3. Finite-Frequency Response 

We now consider the response generated in the disk by 
those terms in the perturbing potential with finite forcing 
frequency. When fl p ^ 0, the second-order differential 
equation for g obtained from equations (|lj), (|l5|), (16) 
and (n8h takes on the general form: 



+ B ~ + <~'<l = S 



'dr 



with 



r(n-np) 



B = fj.fl 



Do 



dln/j, \ 



K 

w 2 



(27) 



(28) 



ft 2 



n _ ftp 

n I ft 2 



-A 



1 d /D x rft 



,02 



5 = i(ft-ft p )X, 



(29) 



(30) 



where we have set D\ = 1 — O p /0 and Z?2 = 
(1 — fip/fi — ?'a v ) 2 — n 2 /fl 2 . By introducing the function 
ft such that: 



1 dft _ 1 / cL4 
ft, dr A \ dr 



equation ( p7| ) can be recast in the form: 

A ( hA*jL) + hC g = hS, 
dr V ar / 



(31) 



(32) 



the solutions of which are given by the Green's function 
method: 



q (r) = — 



9i{r) \ gi{r')s{r')dr' 



+92 (r) / gi (/) s (/) dr' 



(33) 



Here s = h (ft — ftp) frT,fl 2 , is the constant 



„ / i , d 9i dg x 
W=hA[ gi- —32 

ar ar 



(34) 



and 31 and 32 are the solutions of the homogeneous differ- 
ential equation that satisfy the outer and inner boundary 
conditions respectively. In the numerical calculations pre- 
sented below, 31 and 32 are calculated directly from equa- 
tions (|l4|), ( |l5| ) and ( |l8| ) using the fourth-order Runge- 
Kutta procedure with adaptative stepsize control given by 
Press et al. (1986). Computation of the function h then 
allows 3 to be calculated from (|33|). 



5.4. Boundary Conditions 

Equation fl27| ) has a regular singularity at the outer edge 
r = R (see Papaloizou & Terquem 1995). Our outer 
boundary condition is thus that the solution be regular 
there. 

We take the inner boundary to be perfectly reflec- 
tive. The radial velocity in spherical polar coordinates, 

1 /2 

• 2 1 " , thus vanishes at the locations 



( rvi + zv, ) I [r + z 

„2 



/(r 2 



+ zJ 



where 



is the disk inner 



(r, z) such that r 2 
radius. Using the expression (16) of v' z , this means that 
at these locations: 



fig 



1 - ft/ftp 



and then, by continuity: 



fig 



1 - fl/flp 



= 0, 



0, 



(35) 



(36) 



at the inner boundary. 

We note that for the largest values of a v , we expect the 
waves to be damped before reaching the center of the disk, 
such that the results are independent of the inner bound- 
ary condition. For a v = 0, since there is no dissipation at 
the boundaries, we expect no torque to be exerted on the 
disk and no net angular momentum flux to flow through 
its boundaries. 

5.5. Domain of Validity of This Analysis 

In this analysis, we have neglected the variation of 3 with 
z. As noted above, Papaloizou & Lin (1995) have shown 
that this approximation is valid when the scale of variation 
of 3 with radius is larger than the disk semi-thickness H , 
i.e. when 



ldg 

9 dr 



< 



1 

H' 



(37) 



When this is not satisfied, dispersive effects, which have 
been neglected here, should be taken into account. 

Since we also consider linear waves, our analysis is valid 
as long as the perturbed velocities are smaller than the 
sound speed. When they become supersonic, shocks occur 
that damp the waves (Nelson & Papaloizou 1998). How- 
ever, it is not correct to compare the Eulcrian perturbed 
velocities given by the equations (O), ( fl5| ) and ([l6|) with 
the sound speed to know whether the waves are linear or 
not. This is because the system should be invariant un- 
der the addition of a rigid tilt (i.e. a constant) to g. An 
increase of 3 corresponding to the addition of a rigid tilt 
could make the Eulerian perturbed velocities larger than 
the sound speed without the waves getting nonlinear. In 
principle, as pointed out by Papaloizou & Lin (1995), to 
get perturbed quantities that do not depend on the ad- 
dition of a rigid tilt (i.e. that depend only on the radial 
gradient of 3), we should use Lagrangian rather than Eu- 
lcrian variables. However, since the Eulerian perturbed 
radial velocity in spherical polar coordinates depends only 
on the gradient of 3 in the limit of small perturbing fre- 
quency, we can compare this quantity to the sound speed 
to decide whether the motion is subsonic. The criterion 
that has to be satisfied is then: 



8 



TERQUEM 



rv r 



zv' 



(r 2 + z< 



\l/2 



ar 



< C s 



(l-Q p /Q-ia v ) -n 2 /n 2 

(38) 

where we have assumed that £l p <C f2 and z -c r. We 
have kept Q p /Q in the resonant denominator because, even 
though this term is small compared with unity, it may still 
be larger than 1 — n 2 /il 2 . 

The WKB dispersion relation for bending waves gives 
= {c s /2)k where k ~ 1/R is the wavenumber (Pa- 
paloizou & Lin 1995). For hnite perturbing frequency, we 
can then rewrite the above criterion in the form: 



, H 

dg < max — a v 

r 



(39) 



where Sg is the variation of the tilt angle across the disk 
and we have used the fact that 1 — K 2 /fl 2 ~ H 2 /r 2 < H/r. 
By H/r we mean the aspect ratio anywhere in the disk, 
since this quantity does not vary significantly throughout 
the disk. For f2 p = 0, the condition we get is: 



,' H2 

dg < max [ — a v 



(40) 



As pointed out in § 5.2, depending on whether a v is larger 
or smaller than 1 — k 2 /Q 2 , the imaginary part of g is larger 
or smaller than its real part, respectively. Therefore, the 
above criterion and equation (|2l| ) tell us that the variation 
across the disk of the relative vertical displacement pro- 
duced by the secular perturbation is limited by a v along 
the x-axis and by H 2 /r 2 along the y— axis. 
We note that in an inviscid disk, the zero-frequency per- 
turbation produces Sg ~ \lu p \ /tt(R) (Papaloizou, Ko- 
rycansk y fc Terquem 1995). Since, as we have pointed 
out in § 5.2, our analysis is valid only when \lu p \ /fi(ii) is 
smaller than the maximum of H 2 /r 2 and a v , the above 
criterion is always satisfied in an inviscid disk. 

The most likely situation in a protostellar disk corre- 
sponds to a v < H/r. In that case, when the perturbed 
velocities are close to the sound speed, the secular pertur- 
bation produces a steady (in the precessing frame) tilt the 
variation of which across the disk is on the order of H 2 /r 2 
or a v , whichever term is the largest. Superposed on this 
steady tilt, bending waves produced by the finite frequency 
terms propagate, corresponding to a tilt the variation of 
which across the disk is on the order of H/r. This is ex- 
actly what was observed in the SPH simulations performed 
by Larwood & Papaloizou (1997), in which a v was larger 
than H 2 /r 2 . We note that since in the case of protostellar 
disks H/r is likely to be close to 0.1, the variation of the 
vertical displacement across the disk due to the nonsecular 
tidal perturbations can be as large as about a tenth of the 
disk radius while the perturbation remains linear. 



6. ANGULAR MOMENTUM EXCHANGE 

Tidal perturbation of a disk may lead to an exchange of 
angular momentum between the disk and the perturber. 

If the disk is inviscid and does not contain any corotation 
resonance, the nature of the boundaries would determine 
whether such an exchange takes place or not (see Lin & 



Papaloizou 1993 for example). Because of the conserva- 
tion of wave action in an inviscid disk, waves excited at 
the outer boundary propagate through the disk with an 
increasing amplitude if the disk surface density increases 
inward or is uniform (see Lightill 1978 for example). It 
is usually assumed that they become non linear before 
reaching the center and are dissipated through interaction 
with the background flow. Thus, the inner boundary can 
be taken to be dissipative (Papaloizou & Terquem 1995). 
Dissipation at the boundary then introduces a phase lag 
between the perturber and the disk response, enabling a 
net torque to be exerted by the perturber. This torque 
is transferred to the disk through dissipation of the waves 
at the boundary. Because of the conservation of angular 
momentum, the net torque is equal to the difference of an- 
gular momentum flux through the disk boundaries. This 
flux is constant (independent of r) inside the disk, since 
there is not dissipation there. 

The situation is different when a corotation resonance 
is present in the disk, since this singularity provides a 
location where angular momentum can be absorbed or 
emitted (Goldreich & Tremaine 1979; Goldreich & Nichol- 
son 1989). 

When the disk is viscous, its response is not in phase 
with the perturber. A net torque is then exerted on the 
disk even if the boundaries are reflective, and the angular 
momentum flux inside the disk is not constant, since the 
perturbed velocities are viscously dissipated. This situa- 
tion has been investigated by Papaloizou & Pringle (1977) 
and Papaloizou & Lin (1984) in the context of coplanar or- 
bit. Here we present the analysis corresponding to a non 
coplanar orbit. 

We note that whenever the perturber rotates outside the 
disk, the torque exerted on the disk is negative. Through 
dissipation of the waves, the disk then loses angular mo- 
mentum. We first derive an expression of the torque in 
terms of <?, and then relate it to a v using angular momen- 
tum conservation. 

6.1. Expression for the Torque 

The net torque exerted by the perturber on the disk is 
given by: 



T = / Re 
'v 



xRe 



i(tp—n p t) \ 



(41) 



where the integration is over the volume V of the unper- 
turbed disk. Because of the ^-periodicity, the first-order 
term is zero, and the z-component of the torque is then 



T, 



{ Re 




Re 


IV 







d_ 

dip 



(#'< 



i(ip-Qpt) 



which is equivalent to 
1 



T z = ^Im^l^p'**'dV 



dV, 
(42) 

(43) 



Using equations ([12]) and (|l9|), p' in the above equation 
can be replaced in terms of g. The expression of T z then 
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involves 3 ', which, in a Keplerian disk, is equal to ££! 2 . 
We then get 



T z = 7T 



Im (<?) fEr 3 dr 



(44) 



where / is the real amplitude defined by (^) (which has to 
be modified when £l p = in order to take into account the 
disk precession, see § d.z). In the numerical calculations 



presented below, the torque T z will be computed from this 
expression. 

If Jd is the total angular momentum of the disk, the 
tidal effects we have considered would remove the angular 
momentum content of the disk on a timcscalc 



Jd 

T z 



4 My/GM P R 

5 T z ' 



(45) 



We shall compare this timescale with the viscous timescale 
tu,h ~ (t '/ 'i?) 2 r2(i?)~ 1 / 'a h where we will take H/r = 
(H/r) max (for our disk models, H/r does not vary by more 
than 20% in the main body of the disk, as indicated in 
Fig. |l|). Here «h is the viscous parameter that couples to 
the horizontal shear. 

We also define the torque integral T z (r), which is the 
torque exerted between the inner boundary and the ra- 
dius r: 



(46) 



T z (r)=n Im{g)fi:r A dr. 



6.2. Angular Momentum Conservation 

To get the angular momentum conservation equation, 
which is of second order in the perturbed quantities, we 
first take the Lagrangian perturbation of the equation of 
motion (^|) and multiply by the unperturbed mass density. 
To first order, this leads to (Lynden-Bell & Ostriker 1967): 



^ + tf^ + nx(nx«) 



2n 



d 2 $ 

dipdt 



2ttx§ 

at 



dip 



= P A 



-IvpY 

p J 



-F^-pW-p (£ • V) V*. 

(47) 

Here, £ is the Lagrangian displacement vector, A denotes 
the Lagrangian change operator, and we have used the re- 
lation between Lagrangian and Eulerian perturbations of 
a quantity Q to first order in £: 



AQ = Q' + (£ ■ V) 



(48) 



The vector Jl = fife, with k being the unit vector in the 
z-direction. To derive equation (fl7|), we have used the 
fact that at equilibrium F„ = 0. 

We now take the scalar product of equation ( |47| ) with £*, 
which is the complex conjugate of £. Since afot = —i£l p 
and 8/ dip — i, all the terms on the left hand side of the 
resulting equation are real. We then take the imaginary 
part of this equation and integrate it over the volume V 
of the unperturbed disk. This leads to: 



Im {Iv [e ' V{i) + e ' F ' v 



-pC • V*' - p$,* ■ (| • V) V*] dV} = 0, 



(49) 



where we have defined the linear operator V such that 
V (£) = pA (— VP/p) . We note that if the terms on the 
left hand side of equation ( [47| ) had non zero contribution 
to equation J49|), it would be related to the total time 
derivative of the perturbed specific angular momentum of 
the disk. 

Using the perturbed equation of state (|l^) in the form 
AP/P = (1 + 1/n) Ap/p, it can be shown (Lynden-Bell & 
Ostriker 1967) that: 

/ • V (£) dV = B + f C ■ O (£) dV, (50) 
Jv Jv 

where O is a self-adjoint linear operator, and B is a bound- 
ary term: 



B = 



ndS. 



(51) 



The integration is over the surface S of the disk, and n 
is the unit vector perpendicular to this surface, oriented 
toward the disk exterior. Since P and p vanish at r = R 
and z = ±H and we consider models for which r.;„ ~ 0, 
we have B = 0. In addition, since O is self-adjoint, the 
integral on the right hand side of eq uation ( p0[ ) is real. 
The term involving V in equation (B9J) is then zero. 

Using the fact that d 2 ^/drdz = cr 1 i>/dzdr, it can also 
be shown that the term involving V$ is zero. 

Since equation (|4^) has no <p and i-dependence, we can 
now use the tilded variables again. We define £ such that 
£(r, <p, z, t) = £(r, z)exp [i (<p — Q p t)]. We define £* simi- 
larly. The term involving iff' in equation ( [49]) can be writ- 
ten as follows: 



Im 



Im (J p£* ■ ^^' dv 



dVj +Im l^J V'p'*dVj , (52) 

where an integration by parts has been performed and the 
linearized continuity equation: 



p' = -V 



(*-) 



(53) 



(see, for example, Tassoul 1978) has been used. The first 
term on the right-hand side of equation (p3) is zero be- 
cause p vanishes at r = R and z = ±H and ri n ~ 0. 
Equation (M9) then becomes: 



Im ( J p'*WdV 



I in | I £* - F' v dV 



(54) 



From equation (|43]), we see that this is equivalent to: 



T 7 = Im tt 




£* ■ F' v rdrdz 



(55) 



H 



We now use the relation between the Eulerian velocity 
perturbation and the Lagrangian displacement: 
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v = § + (v-V)| (56) 



together with the expression (|13| ) of P^, the relation 
and the fact that p (±H) = to transform equation 
into: 



Trim 



pa v r- 



pa v r- 



(v'*/z) 



\v>/z\ 2 + K/z\ 2 



i - n p /n 
%/*) 



dr 



(1 - fip/rt) 



K 

2^2 



(57) 



In the numerical calculations presented below, the 
torque will be computed from equation (Q) and compared 
with the result of the above equation 

In a Kcplerian disk, k 
a f « 1, equation (fLa) is well approximated by 



f2. Then, when £l p <C f2 and 



z 



2il 2 



so that T z can be approximated by: 



3tt 







K 


/ pa v r 






z 



dr. 



(58) 



(59) 



As expected, the torque is negative. Here we consider 
linear perturbations such that the perturbed velocities re- 
main smaller than the sound speed. We have pointed out 
in § 5.5 that, since the system should be invariant under 
the addition of a rigid tilt to g, it is not in principle correct 
to compare the Eulcrian perturbed velocities to the sound 
speed to know whether the waves are linear or not. How- 
ever, since the dominant term in the expression ( |l4| ) of v' r 
is the term proportional to the gradient of g, we can com- 
pare this velocity to the sound speed to decide whether 
the motion is subsonic. Therefore we can get an estimate 
of the maximum value of the torque by setting \v' r \ ~ c s , 
/Lt ~ pH 3 and dropping the integration. This leads to a 
minimum value of t& which is equal to t u ,h X ah/« v - Only 
if a v were larger than ah could td be smaller than t v ,h- 

An other useful quantity for interpreting our results is 
the flux of vertical angular momentum in the radial direc- 
tion at radius r: 



J(r) 



H 



I z=-H Jtp=0 

which can also be written 

J(r) = npr 2 Re 



pr 2 Re {v' v ) Re (v' r ) dipdz, (60) 



_£_ 

Z 



(61) 



The torque T z (r) exerted between the radii Ti n and r re- 
sults in a change of angular momentum of the disk. Part 
of this angular momentum is advected through the bound- 
aries (this is J(r) — J(rj n ) ~ J(r) since Ti n ~ 0), and 
part is dissipated between these radii. Thus the quantity 
T z (r) — J(r) is the amount of angular momentum deposited 
in the disk between the inner boundary and the radius r. 



As we shall see, this quantity is negative, which means 
that the disk loses angular momentum. 

In the numerical calculations presented below, the 
torque integral will be computed from equation (E6I) and 
J(r) will be computed directly from the velocities using 
equations ( |l4| ) and (|Tq). 

7. NUMERICAL RESULTS 

We compute the torque exerted on protostellar disks by 
a perturber with an orbit that is inclined to the plane of 
the disk. For a coplanar system, one expects the circum- 
primary disk to be truncated by tidal effects in such a 
way that its radius is not greater than about one-third of 
the separation of the system (Papaloizou & Pringle 1977). 
Larwood et al. (1996) have shown that tidal truncation is 
only marginally affected by the lack of coplanarity. Thus 
in our models we adopt D > 3i?. Since we consider binary 
mass ratio of unity [M s — M p ), the perturbing potential 
is small compared with the potential of the central star, 
justifying the assumption of small perturbation. 

For computational convenience, we normalize the units 
such that M p = 1, R = 1, and Qk{R) — 1- Apart for 
the calculation of the epicyclic frequency, the rotation law 
in the disk is taken to be Keplerian (the departure from 
the Keplerian law, which occurs mainly at the disk outer 
edge, does not exceed ~ 3%). The inner boundary of the 
disk is ri n = 10~ 4 , and the outer boundary is taken to be 
font = 0.99 rather than 1 in order to avoid a zero surface 
mass density at the outer edge. 

We present here the results of some disk response cal- 
culations for an inclination S = 7r/4, a polytropic index 
n = 1.5 and a disk mass M = 10~ 2 . Results for an arbi- 
trary inclination S < n/2 (prograde orbit) or tt/2 < 5 < ir 
(retrograde orbit) can be obtained straightforwardly since 
T z (and then 1/trf) corresponding to O p = ±2w and is 
proportional to (1 ± cos<5) 2 sin 2 5 and sin 2 26, respectively. 
As shown by equation (24), the precession frequency scales 
with cos 6. Since g is independent of M (see equations jyj, 
pi, @ and @), expression (f3[) shows that T z oc M. 
The quantities u> p , td and t v ^ are independent of M. 

7.1. Uniform a v 

In table |l| we summarize the results obtained for dif- 
ferent values of the parameters, which are the separa- 
tion of the system D, the viscous parameter a v and the 
maximum value of the relative semi-thickness of the disk 
(H/r) max . The quantities we compute are the torque T z 
associated with each of the perturbing frequencies 2lo (pro- 
grade term), — 2u> (retrograde term) and 0, the disk pre- 
cession frequency oj p and the ratio of the tidal timescale td 
to the viscous timescale t v ^. We have taken here ah = a v - 
The two values of the torque have been computed using 
equations (Q) and (|57j), the result corresponding to equa- 
tion (57) being in parentheses. The timescale td is com- 



puted Trom equation 
contributions from eac 



with T z being the sum of the 
i frequency. 



We note that when either one of the criterions (|37|) 

or 
is a 



is not satisfied (cases indicated by a footnote), it 
,Iways close to the disk outer edge. This is because the 
density varies more rapidly there than in the other parts 
of the disk. Also the scale of variation of g never becomes 
smaller than H by less than a factor of a few. Similarly, 
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when the radial velocity in spherical polar coordinates be- 
comes larger than the sound speed, it is only by a factor of 
a few. We note that the results we present in this case can 
still be used since the perturbed quantities can be reduced 
by choosing appropriate scaling parameters, like the initial 
inclination 6 or the mass of the disk M. 

In all cases the con dition |w p | /^(i?) < max {H 2 /r 2 , a v ) 
is satisfied (see § 5^2). The condition for near rigid body 



precession, H/R » \oj p \ /Q(R), is also always satisfied. 

We note that there is a good agreement between the 
values of the torque given by equations (Q) and ([37]) . In 
addition, in the cases we computed, we found this latter 
equation to be very well approximated by ( |59| ) . 

To simplify the discussion, we note T+, T~ and T° the 
torque associated respectively with f2 p = 2lu, —2uj and 
0. Furthermore, we separate the ^-dependence by writing 
Tf = T±(l ± cos 2 sin 2 6 and T° = T sin 2 25 with T± 
and Tq being independent of S. 

Since part of this section is devoted to the tilt g, we re- 
call that g is related to the relative vertical displacement 
through equation (pl|). When Q p — 0, Re(g) and Im(g) 
are the relative vertical displacement along the y and x- 
axis respectively. 

7.1.1. Finite-Frequency Response 

Figure || shows the real and imaginary parts of g versus 
r for models 2 and Q p = 2lu (figures for £l p = — 2lu are 
similar). We see that the wave-like character of the disk 
response disappears when a v becomes larger than H/r 
(which is 0.05 in these models). This is in agreement with 
Papaloizou & Pringle (1983) and Papaloizou & Lin (1995) 
(and also with the relativistic generalization of Demian- 
ski & Ivanov 1997 and Ivanov & Illarionov 1997) who 
have shown that, in a near Keplerian disk (self-gravitating 
or not), the longest wavelength disturbances undergo a 
transition between wave-like and diffusive behavior when 
a v ~ H/r. We note, however, that when a v > H/r, there 
are still some oscillations in the half outer part of the disk. 
This is because bending waves have a long wavelength, 
such that they can still penetrate relatively far before be- 
ing diffused out. 

The magnitude of Re{g) does not vary with a v , since it 
is controlled mainly by the radial pressure force. On the 
contrary, the magnitude of Im(g) is much more sensitive 
to the viscosity, although it does not vary monotonically 
with a v . 

By comparing models 2a and 3c and also models 2c and 
4c, which differ only by the value of {H/r) max , we have 
checked that the wavelength of the response is proportional 
to H/r. This is expected since bending waves propagate 
with a velocity which, being half the sound speed, is pro- 
portional to H/r, and the wavelength is proportional to 
the wave velocity. 

In all the cases we have computed, except model 3f, 
T+ > T~ , but T_ > T + , in agreement with Papaloizou & 
Terquem (1995). In general, the shorter the wavelength 
of the response, the smaller the torque, and the wave cor- 
responding to the retrograde term has a slightly longer 
wavelength than that corresponding to the prograde term. 
The difference between T + and 7L (and To) is still sig- 
nificant for the largest values of the separation computed 
here. We expect all these quantities to converge toward 



the same value when 2u> approaches 0, which probably 
means here 2u)/£l(R) <^ a v since the torque is controlled 
by viscosity. This is not satisfied even for model 3i which 
has 2u)/Cl(R) = 0.09. However, we observe that the differ- 
ence between these torques is reduced from a factor 9 to a 
factor less than 4 when model 3i is run with a v — 0.1. 

In Figure || the net torque T z is plotted versus D in 
a semi-log representation for a v = 10~ 3 (models 3 plus 
other values of D) and a v = 10~ 2 (models 4 plus other val- 
ues of D) and for fl p = ±2uj and 0. Before commenting on 
the resonances, we note that the torque associated with the 
finite frequency terms out of resonance does not decrease 
when the separation increases up to O ~ 8. This con- 
firms the results of Papaloizou & Terquem (1995). There 
are competing effects acting when the separation of the 
system is increased. On the one hand, the wavelength of 
the tidal waves becomes larger (since it is inversely pro- 
portional to its frequency), which tends to increase the 
torque. On the other hand, the amplitude of the perturb- 
ing potential decreases, which tends to reduce the torque. 
It seems that the first effect is dominant for D < 8. It 
seems like the torque begins to decrease with the separa- 
tion when the wavelength becomes comparable to the disk 
radius, which happens for D ~ 9. This is reasonable since 
the wavelength cannot increase further. 

In Figure @ we can see resonances, which occur when 
the frequency of the tidal waves is equal to that of some 
free normal global bending mode of the disk, and cause 
the torque to become very large, even infinite if there is 
no dissipation (Papaloizou & Lin 1984). In a resonance, 
the torque indeed increases when a v decreases, in contrast 
to what is observed out of resonance. It seems that there 
is no longer any resonance for D > 10. This would not 
be surprising since the wavelength of the response is then 
comparable to the disk radius. A more detailed descrip- 
tion of the resonance which occurs at D ~ 8.5 is given in 
the Appendix. 

Since the resonances depend on the spectrum of the 
free normal bending modes of the disk, they are model- 
dependent. In particular, they occur only if the the waves 
can be at least partially reflected at the inner boundary, 
otherwise there is no free normal global bending mode in 
the disk (see also Hunter & Toomre 1969 for the case of a 
purely self-gravitating disk). 

The values of the torque displayed in Figure || corre- 
spond to 5 — 7r/4. However, for this particular value of 
6, the perturbed velocities may become larger than the 
sound speed close to a resonance. These values of the 
torque should then be scaled to be used for other values 
of 5, or for other disk parameters. For a v = 10~ 2 , the 
resonances are very weak, indicating that the waves are 
almost completely dissipated before reaching the disk in- 
ner boundary. 

We have compared our results with those obtained by 
Papaloizou & Terquem (1995) for an inviscid disk in which 
the inner boundary is dissipative. Remembering that they 
had 6 — 7r/2, we found that, in general, we could repro- 
duce their results with a v ~ 10~ 2 . This confirms that, 
when a v ~ 10~ 2 , the waves are dissipated before reaching 
the disk inner boundary. 

Out of resonance, the torque Tf- associated with the fi- 
nite frequency terms is proportional to a v as long as this 
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parameter is smaller than some critical value. Then, for 
the larger values of a v , Tf- is independent of the viscosity. 
This is in agreement with the theoretical expectation (Gol- 
dreich & Tremaine 1982) and the work of Papaloizou & 
Terquem (1995). We indeed expect the torque to be inde- 
pendent of whatever dissipation acts in the disk providing 
the waves are dissipated before reaching the disk inner 
boundary. If the waves are reflected on the disk inner 
boundary, the situation is more complicated because the 
angular momentum carried by the waves has a different 
sign depending on whether the waves are in-going or out- 
going. Therefore part of the angular momentum lost by 
the disk while the waves propagate inwards is given back 
by the waves propagating outwards. The net amount of 
angular momentum lost by the disk then depends on how 
efficiently the waves are reflected. Using this argument, 
we then find that, when H/r = 0.1 and 0.05, the waves 
are damped before reaching the disk inner boundary when 
a v is larger than ~ 5 x 10~ 2 and ~ 5 x 10 -3 , respectively 
(for H/r = 0.1, this is a bit larger than the value found 
above) . Such a dependence of this critical value of a v with 
H/r is expected. Indeed, the shorter the wavelength, the 
more easily are the waves dissipated. Since the wavelength 
is proportional to H/r (see above), we then expect that 
for a given a v dissipation is more efficient for small values 
of H/r. For the same reason, the net torque Tf- increases 
with H/r (see models 3c and 2a and models 4c and 2c). 
The coupling between the response and the perturbing po- 
tential is indeed more efficient when the response has a 
long wavelength, and the integral in equation (Q) from 
which the torque is calculated is then larger. 

7.1.2. Zero -Frequency Response 

Figure |] shows the real and imaginary parts of g versus 
r for models 2 and Q, p = 0. As expected, we observe that 
the zero-frequency response is not wave-like. 

From equation (|2^) we see that, as long as a 2 <C 
l-K 2 /ft 2 ~ H 2 /r 2 , which is the case for a v < 10 -2 , Re(g) 
is almost independent of the viscosity. When a v = 0.1, a 2 
dominates over 1 — K 2 /f2 2 , and since it appears in the ex- 
pression of Re(g) with the minus sign, Re(g) becomes neg- 
ative. The particular behavior of Re(g) in the outer parts 
of the disk is produced by the fact that the surface mass 
density drops to zero there, thus increasing 1 — K 2 /fl 2 . As 
expected from equation (|2q), Im(g) is proportional to a v . 

Papaloizou, Korycansky & Terquem (1995) computed 
the real part of g for an inviscid disk with an equilibrium 
state similar to that we have set up here and for fl p = 
and D = 7. We have checked that the results we get for 
the smallest values of a v are in complete agreement with 
theirs. 

Equation ( |44] ) shows that T z depends only on Im(g), not 
on Re(g). Since for the secular response Im(g) oc a v , T° 
is also proportional to a v . This is borne out by the numer- 
ical results. Also equation ( p6| ) indicates that g oc l/H 2 . 
Thus if we vary H/r while keeping S constant, we get 
from equation ( p4| ) that T z oc l/H 2 . This is confirmed 
by models 3c and 2a and models 4c and 2c, which differ 
only by the value of (H/r) max . This is in contrast to the 
finite frequency response, for which the torque increases 
with H/r (see above). 

Figure || indicates that T° 3> Tf- out of resonance. Since 



we also have Tq T± , T~ is always going to be much larger 

, o 

than T z for 5 not too close to 7r/2 (typically for 5 < 70 

or 110° < 8 < 180°). 

In contrast to the finite frequency response, the zero- 
frequency response is not wave-like, and we then expect 
T° to decrease continuously with D. This is indeed what 
we observe in Figure pi 

7.1.3. Angular Momentum Dissipation as a Function of 

Radius 

Figures ^ and ^| show T z (r) — J(r) normalized to unity 
versus r for various models. This represents the angu- 
lar momentum deposited in the disk between the inner 
boundary and the radius r. The fact that this quantity is 
negative means that the disk loses angular momentum. Of 
course, \T z (r) — J(r)\ should increase with r. However, we 
observe in Figure H that for D = 3 and fi = 2lo this func- 
tion has a little 'hump' close to the disk outer edge. We 
interpret it as being a numerical effect, since its amplitude 
is found to depend on the accuracy of the calculations. 

We first observe that bending waves are able to trans- 
port a significant fraction of the negative angular momen- 
tum they carry deep into the inner parts of the disk. For 
the parameters corresponding to model lb, Figure]^ shows 
that the retrograde and prograde waves transport respec- 
tively 35% and 60% of the total negative angular momen- 
tum they deposit into the disk at radii smaller than 0.3 and 
0.4 respectively. In contrast, only 20% of the total negative 
angular momentum carried by the zero-frequency pertur- 
bation is deposited at radii smaller than 0.4. We note 
however that, out of resonance, the net amount of angular 
momentum deposited at small radii by the fl p — per- 
turbation is not necessarily smaller than that deposited 
by the finite frequency perturbations. Also, the fraction 
of angular momentum deposited close to the disk outer 
edge by the finite frequency responses (especially the ret- 
rograde wave) is larger than that deposited by the secular 
response. 

It appears that the secular response gets dissipated very 
progressively throughout the disk, whereas the prograde 
and retrograde waves are dissipated preferentially at some 
locations in the disk. The position of these locations seems 
to depend on the wavelength of these waves. 

The resonances do not affect the distribution of dissi- 
pation of angular momentum in the disk, although they 
obviously increase the absolute value of the actual amount 
of angular momentum deposited. This is due to the fact 
that the magnitude of the response, not its wavelength, is 
affected by a resonance. 

In Figure ^ we have plotted T z (r) — J(r) versus r for 
Clp = —2oj and different values of a v . We see clearly on 
these plots the transition that occurs when a v ~ H/r. 
When a v becomes larger than H/r, the wave-like response 
of the disk disappears. Thus, the curves corresponding to 
Q p = —2co have a shape similar to those corresponding to 
the secular response. 

In general, when a v is increased, the fraction of nega- 
tive angular momentum deposited at small radii decreases. 
This is expected since the waves, excited predominantly in 
the outer parts of the disk, get dissipated more easily at 
large radii when a v is large. 
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We observe in Figure |6[a that there is hardly any differ- 
ence between the cases a v = 1(T 4 and a v = 1CT 3 . This 
means that whatever the amount of angular momentum 
transported by the wave, the same fraction is deposited be- 
tween ri n and r in both cases. This probably indicates that 
the wave is able to reach the disk inner boundary for these 
values of a v . For larger a v for which the wave is completely 
dissipated before reaching the inner boundary, no angular 
momentum is deposited at small radii, so that the curves 
look different. This indicates that, for (H/r) max =0.1, 
the waves do not reach the inner boundary when a v is 
larger than - 1CT 2 . For (H/r) max = 0.05 (Figure |.b), 
the critical value of a v above which the waves do not reach 
the inner boundary is between 10 -3 and 5 x 10~ 3 . These 
values arc in agreement with those found above. Our cal- 
culations also show that a larger fraction of the negative 
angular momentum is deposited at small radii when H/r 
is decreased. However, we have checked that this effect 
is less significant when the distance interior to the outer 
boundary occupied by the taper of the surface mass den- 
sity is reduced along with the semi-thickness. 

7.2. Nonuniform a v 

We have run model lc with a nonuniform viscosity 
a v = 10 _2 a; 5 and a v = 10~ 2 (1 — x) 5 . These models cor- 
respond to a disk where respectively the inner and outer 
parts are nonviscous. 

Figure ^.a shows the imaginary part of g versus r for 
O p = corresponding to these viscosity laws. The uni- 
form a v case is also plotted for comparison. The real part 
of g is not shown since it does not depend on a v for these 
small values of the viscosity (see discussion above). We 
see that when the outer parts of the disk are nonviscous, 
Im(g) in the inner parts is similar to the uniform a v case, 
while it has an almost constant value in the outer parts. 
When the inner parts of the disk are nonviscous, Im(g) 
is almost zero there, while it is significant in the outer 
parts. This behavior can be understood by remembering 
that Im[g(r)] is the integral from ri n to r of a function 
proportional to a v (see equation pq| ). Since Im(g) varies 
less through the disk than in the uniform a v case, v' r given 
by equation (|l4| ) is smaller. This, together with a smaller 
a v , leads to a smaller net torque given by equation (|59|). 

When Sip ^ 0, there is not such a difference between the 
uniform and nonuniform a v cases. The torque is smaller 
in the nonuniform case, but Im{g) does not change signif- 
icantly from one case to the other. The reason is probably 
that Im(g) does not vary monotonically with a v , as shown 
in Figure b[ 

Figure 0Tb shows T z {r) — J(r) versus r for the same mod- 
els and fl p = 0. We observe that, as expected, no angular 
momentum is deposited at radii where a v is very small 
(small radii for a v oc x 5 and large radii for a v oc (1 — a;) 5 ). 

8. DISCUSSION 

In this paper, we have calculated the response of a vis- 
cous gaseous disk to tidal perturbations with azimuthal 
mode number m = 1 and odd symmetry in z with respect 
to the equatorial plane for both zero and finite perturbing 
frequencies. We concentrated on these types of perturba- 
tions because they arise for inclined disks and they lead to 



a long-wavelength response. The effects of perturbations 
with even symmetry, which occur in the coplanar case, 
may be superposed on the effects of those with the odd 
symmetry. 

Since the response of a viscous disk is not in phase with 
the perturbing potential, a tidal torque is exerted on the 
disk. When the perturber rotates outside the disk, this 
torque results in a decrease of the disk angular momen- 
tum. Part of the (negative) angular momentum of the 
perturbation is carried by tidal waves away from the lo- 
cation where the torque is exerted, and part is dissipated 
locally through viscous interaction of the waves with the 
background flow. 

We have shown that the tidal torque is comparable to 
the horizontal viscous stress acting on the background flow 
when the perturbed velocities in the disk are on the order 
of the sound speed c s . If these velocities remain subsonic, 
the tidal torque can exceed the horizontal viscous stress 
only if the parameter a v which couples to the vertical 
shear is larger than the parameter ah, which is coupled 
to the horizontal shear. We note that, so far, there is no 
indication about whether these two parameters should be 
the same or not. Nelson & Papaloizou (1998) have found 
that, when the perturbed velocities exceed c s , shocks re- 
duce the amplitude of the perturbation such that the disk 
moves back to a state where these velocities are smaller 
than c s . When shocks occur, the tidal torque exerted on 
the disk may become larger than the horizontal viscous 
stress. 

We have found that, in protostellar disks, bending waves 
are able to propagate and transport a significant fraction of 
the negative angular momentum they carry in the disk in- 
ner parts. This is due to their relatively large wavelength. 
Therefore, tidal interactions in noncoplanar systems may 
not be confined to the regions close to the disk outer edge 
where the waves are excited. For the disk models we have 
set up, which extends down to the stellar surface, the value 
of a v above which the waves get dissipated before reaching 
the disk inner boundary varies between 5 x 10~ 3 and 10~ 2 
for a disk aspect ratio between 0.05 and 0.1 (this critical 
value of a v increases with the disk semi-thickness). We 
note that if we had set up an annulus rather than a disk, 
the surface mass density would drop at the inner bound- 
ary. Then, for small a Vl the waves would have a tendency 
to become nonlinear and then to dissipate before reaching 
the inner boundary. Thus, the behavior of bending waves 
could be similar in an annulus with small a v and in a disk 
with larger a v . In the limit of small viscosity, dissipation 
could also occur as a result of parametric instabilities (Pa- 
paloizou & Terquem 1995; C.F. Gammie, J.J. Goodman & 
G.I. Ogilvie 1998, private communication), the effect of 
which may be to lead to a larger effective a v . However, 
since the waves we consider here propagate in a nonuni- 
form medium, it is not clear whether they would be effi- 
ciently damped by these instabilities. 

It the disk model allows at least partial reflection from 
the center, the tidal interaction becomes resonant when 
the frequency of the tidal waves is equal to that of some 
free normal global bending mode of the disk. For the equi- 
librium disk models we have considered, we have found 
that a particularly strong resonance occurs when the sep- 
aration is about 8.5 times the disk radius. The torque 
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associated with the finite frequency terms then increases 
by many orders of magnitude. However the response may 
be limited by shocks and nonlinear effects. The strength of 
the resonance is inversely proportional to a v . In our calcu- 
lations, the range of separation for which the torque is sig- 
nificantly increased is rather large, which means that the 
effects of tidal interaction, like disk truncation, may occur 
even when the separation of the binary is large. In addition 
to this strong resonance we have found a few other reso- 
nances with different strengths and widths. From these 
calculations, it appears that the probability for the inter- 
action to be resonant may be significant. However, these 
results depend on the particular equilibrium disk models 
we have set up. 

We note that the disk is expected to be truncated such 
that the inner Lindblad (2:1) resonance with the compan- 
ion, which could provide effective wave excitation, is more 
or less excluded from the disk. 

Out of resonance, we find that the torque associated 
with the zero-frequency perturbing term is much larger 
than that associated with the finite frequency terms. Of 
course, if the separation of the system is large enough 
(at least 10 times the disk radius in our models), the dif- 
ference between finite frequency and secular terms disap- 
pears. When the secular response is dominant, the tidal 
torque is proportional to a v in the limit a T <l. This has 
the consequence that the ratio of the tidal timescale td 
(time that would be required for the tidal effects we have 
considered to remove the angular momentum content of 
the disk) to the disk viscous timescale t v ^ is proportional 
to ah/a v . 

What observable effects would these tidal interactions 
produce in protostellar disks ? First, as we mentioned 
in the introduction, they would lead to the precession 
of the jets that originate from these disks. The preces- 
sion timescale that we can infer from the observation of 
jets that are modeled as precessing is consistent with this 
motion being induced by tidal effects in binary systems 
(Terquem et al. 1998). 



Also, as we have already pointed out in § 5.5, the sec- 
ular perturbation produces a tilt the variation of which 
can be up to a v along the line of nodes and H 2 /r 2 along 
the direction perpendicular to the line of nodes. Super- 
posed on this steady (in the precessing frame) tilt, there 
is a tilt produced by the finite frequency perturbations, 
the variation of which across the disk can be up to H/r. 
Such asymmetries in the outer parts of protostellar disks 
could be observed. Because the variations along the line of 
nodes and along the perpendicular to the line of nodes de- 
pend on a v and H/r, respectively, observations of warped 
protostellar disks have the potential to give important in- 
formation about the physics of these disks. 

Protostellar disks are believed to be rather thick, i.e. 
H/r ~ 0.1. Our calculations show that for such an aspect 
ratio, the viscosity above which bending waves are damped 
before reaching the central star is a v ~ 10~ 2 . Observa- 
tions seem to indicate that ah in protostellar disks is on 
the order of 10 _3 -10 -2 . If a v is smaller than or compa- 
rable to 10 -2 , resonances may occur, as described above, 
providing the disk inner edge allows some reflection of the 
waves. In that case, we could observe truncated disks with 
sharp edges even when the binary separation is large. If 



a v is larger than but close to 10~ 2 , bending waves are still 
able to propagate throughout most of the disk. In addi- 
tion, given that protostellar disks are rather thick, these 
waves propagate fast (they cross the disk on a time com- 
parable to the sound crossing time). In that case, it may 
be possible to observe some time-dependent phenomena 
with a frequency equal to twice the orbital frequency. 

In the case of protostellar disks, comparison between 
observations and theory is now becoming possible. A non- 
coplanar binary system, HK Tau, has been observed for 
the first time very recently (see § 1.1), and subsequent 



work will be devoted to interpreting these observations. 
In the meantime, we can comment briefly on some obser- 
vations related to X-ray binaries. 

There is evidence from the light curve of X-ray binaries, 
such as Hercules X-l and SS 433, that their associated ac- 
cretion disks may be in a state of precession in the tidal 
field of the binary companion. Katz (1980a, 1980b) has in- 
dicated that the observed precession periods of these two 
systems are consistent with the precession being induced 
by the tidal field of the secondary. 

In the case of SS 433, it is interesting to note that an 
additional "nodding" motion with a period half the or- 
bital period is observed (see Margon 1984 and references 
therein). Katz et al. (1982) suggested that this nodding 
motion is produced by the gravitational torque exerted by 
the companion. In their model, the nodding is viscously 
transmitted from the outside, where the torque is signifi- 
cant, to the interior, where the jets are seen to respond to 
the motion. This implies an extremely large viscosity in 
the disk. 

Here we point out that the observed period of this mo- 
tion is consistent with the nodding being produced by 
bending waves. As these waves propagate with the sound 
speed and the disk is observed to be very thick (Mar- 
gon 1984), they could transmit the motion through the 
disk on a timescale comparable to its dynamical timescale. 
For the waves to reach the disk interior, a v would have to 
be smaller than H/r. It is not clear that this condition 
is satisfied here. The disk viscosity would then probably 
have to be smaller than that predicted by the model of 
Katz et al. (1982). 

We finally observe that a rapid communication time 
through the disk also makes problems with the radiation- 
driven warping model proposed by Maloney & Begel- 
man (1997) to explain the precession of the disk in SS 433. 
In this model, the communication between the different 
parts of the disk, which is assumed to occur because of 
viscosity alone, is on a timescale that is characteristic for 
mass to flow through the disk. It is therefore very slow 
since the viscous timescale must be long for the warping 
instability to occur (Pringle 1996). 
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APPENDIX 

We describe here in more detail the resonance that occurs at D ~ 8.5, since it is particularly strong. For fi p = 2u> and 
— 2u>, resonance occurs at D — 8.25 and D — 8.425 respectively. We observe that, for the values of a v considered here, 
the range of D for which the torque is significantly increased is rather large. We see that the strength of the resonance is 
proportional to a^ 1 - This is expected since the strength of the resonance is inversely proportional to the damping factor, 
which itself is proportional to a v when dissipation is due to viscosity. We also note that the integral of the torque over the 
range of frequencies in this resonance is independent of a v . This can be understood in the following way. For Q p = 2lo 
(the same argument would apply for f2 p = —2a;), the torque at or close to resonance can be written 

_ -Aa v 

2 — 2 ' 

4 (u - u> ) + al 

where ui — uj/Q(R), uj is the resonant frequency and A is an amplitude which does not depend on a v . We have calculated 
that A ~ 8 x 1CT 11 . For the values of a v considered, we have checked that the integral of the torque from d>o — Aw to 
ujq + Ad) does not depend very much on whether 2Ad> is taken to be the range of frequencies in the resonance or for which 
expression (|l|) is valid (the latter being smaller). We then approximate the integral I of the torque over the resonance by 
using (Q), so that 

I = -A Arctan^^-. (2) 

For a v = 1CP 3 or 1CP 4 , the range of frequencies 2Ad) for which (|l|) is valid is large compared with a v , so that / ~ —An/2, 
independent of a v . For a v = 10~ 2 , although 2Ad> is only a few times a v , this relation is still valid within a factor less 
than 2 (due to the rapid convergence of the function Arctan) . 
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Note. — Listed are the separation of the system D, the viscous parameter a v , the maximum value of H/r, the torque T z associated with each 
of the perturbed frequencies 2u>, —2io and 0, the disk precession frequency lu p and the ratio of the tidal timescale t a to the viscous timescalc 

"Cases where, at the outer edge of the disk, g varies on a scale smaller than H. 

'Cases where, at the outer edge of the disk, the radial spherical velocity is larger than the sound speed. 

c Cases where, at the outer edge of the disk, g varies on a scale smaller than H and the radial spherical velocity is larger than the sound speed. 
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Fig. 1. — Surface mass density (top left), relative semithickness (top right), ratio of the sound speed to the Kcplerian 
frequency (bottom left) and ratio of the epicyclic frequency to the Keplerian frequency (bottom right) vs. r for the 
equilibrium disk models. The parameters are n — 1.5, M = 10~ 2 M P and (H/r) max = 0.1. 
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Fig. 2. — Real (panel a) and imaginary (panel b) parts of g vs. r for models 2 and fl p — 2lu. The different curves 
correspond to a v = 0.1 (solid line), 5 x 10~ 2 (dotted line), 10~ 2 (short dashed line), 5 x 10~ 3 (long dashed line) and 
10~ 3 (dot-short dashed line). The parameters are D = 6 and (H/r) max = 0.05. The wavelike character of the response 
disappears when a v > H/r. 
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Fig. 3. — Net torque \T Z \ vs. D in a semilog representation for (H/r) max = 0.1. The panels a and b correspond to 
a v = 10~ 3 (models lb and 3) and a v = 1CT 2 (models lc and 4) respectively. The different curves correspond to Q p = 2ui 
(solid line), —2uj (dotted line) and (dashed line). These plots show several resonances. 
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Fig. 4. — Real (panel a) and imaginary (panel b) parts of g vs. r for models 2 and Q p — 0. The different curves correspond 
to a v = 0.1 (solid line), 5 x 10~ 2 (dotted line), 10 -2 (short dashed line), 5 x 10~ 3 (long dashed line) and 10 -3 (dot-short 
dashed line). The parameters are D = 6 and (H/r) max = 0.05. The zero-frequency response is not wave-like. 
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Fig. 5. — T z (r) — J(r) (angular momentum dissipated between n n and r) normalized to unity vs. r for a = 10 3 , D = 3 
(model lb) and Q p = (solid line), 2uo (dotted line) and — 2ui (dashed line). 
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Fig. 6. — T z (r) — J(r) (angular momentum dissipated between r in and r) normalized to unity vs. r for £l p = —2ui. The 
panel (a) corresponds to models 1, i.e. D = 3, {H/r) max — 0.1 and a = 10~ 4 (solid line), 1CP 3 (dotted line), 10~ 2 (short 
dashed line) and 0.1 (long dashed line). The panel (b) corresponds to models 2, i.e. D = 6, (H/r) max = 0.05 and a = 10~ 3 
(solid line), 5 x 10~ 3 (dotted line), 10~ 2 (short dashed line), 5 x 10~ 2 (long dashed line) and 0.1 (dot-short dashed line). 



ACCRETION DISKS AND BENDING WAVES 



23 





Fig. 7. — Imaginary part of g {panel a) and T z (r) — J(r) (angular momentum dissipated between r in and r) normalized 
to unity (panel b) vs. r for fl p — and a nonuniform a. The different curves correspond to a = 1CP 2 (solid line), 
10~ 2 (1 — x)° (dotted line) and 10~ 2 a; 5 (dashed line). The parameters are D = 3 and (H/r) max = 0.1. 



